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SUMMARY 


The laminar flow of a slightly viscous incompressible fluid that 
issues from a slit and passes over a flat plate is investigated in a 
region far enough from the slit for the boundary— layer equations to be 
valid. By assuming similar velocity distributions along the plate, the 
partial differential equation for the boundary layer is reduced to a 
third-order nonlinear ordinary differential equation. This equation is 
integrated by numerical means for the required boundary conditions. The 
solution gives the velocities at points in the fluid and the surface 
friction at points on the plate. 

Some of the specific results obtained are that the velocities paral- 
lel to the plate vary inversely as the 1/2 power of the distance from 
the slit, that the velocities perpendicular to the plate vary inversely 
as the 3/4 power of the distance from the slit, that the width of the 
disturbed region increases as the 3 /h- power of the distance from the 
slit, that the surface friction is independent of the viscosity of the 
fluid and varies inversely as the 5/4 power of the distance from the 
slit, that the rate at which momentum parallel to. the plate passes 
through a^ plane normal to the plate varies inversely as the 1/4 power 
of the distance from the slit to the normal plane, and that the quantity 
of fluid passing through a plane normal to the plate per unit time varies 
directly as the l/4 power of the distance from the slit to the normal plane. 


INTRODUCTION 


The laminar flow of an infinite uniform stream of incompressible 
fluid of small viscosity over a flat plate at zero angle of attack was 
first investigated analytically by Prandtl (reference 1). By a proper 
choice of variables the partial differential equation of the Pr an dtl 
boundary layer was reduced to a nonlinear ordinary differential equation 
of the third order. An accurate solution of the equation that gave the 
skin-friction coefficient for the plate and the velocities in the field 
of flow was first obtained by Blasius (reference 2) and later by others. 
(For example, see references 3 to 6.) The numerical values obtained by 
the different investigators, all of whom used either series or numerical 
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methods to integrate the nonlinear differential equation, were about the 
same. 


The two-dimensional laminar flow of a slightly viscous incompressible 
fluid from a slit into an infinite region of still fluid was investigated 
analytically by Schlichting (reference 7 ) and by Bickley (reference 8 ). 

By introducing new variables into the Prandtl boundary-layer equation 
•and by using the fact that the rate of flow of momentum throug h a cross 
section of the Jet is independent of the distance from the slit, the 
partial differential equation of the boundary layer was reduced to a 
third-order nonlinear ordinary differential equation. Schlichting inte- 
grated the equation by the method of series and Bickley later integrated 
the equation in a closed form. The solutions gave the velocities at 
points in the field of flow.. In addition to the solutions for the flow 
over a flat plate and for the flow from a slit, reference 9 also gives 
the few other known solutions of the boundary— layer equations. The solu- 
tions are given for the flows near the stagnation points of plane bodies 
and of bodies of revolution, for the flow along a wall in a converging 
channel, and for the flow in a round Jet from which the fluid issues 
from a small hole in a wall. Any additional solution of the boundary- 
layer equations Is then, aside from any practical application, of some 
importance in boundary— layer theory. 

The purpose of the present work is to investigate by means of the 
Prandtl boundary— layer equations the laminar flow of a slightly viscous 
incompressible fluid that issues from a slit and passes over a flat 
plate. The arrangement (fig. 1) may also be interpreted as the laminar 
flow into a very large container of still fluid from a slit at the 
intersection of two of the walls,. By use of the Prandtl boundary- layer 
equations together with the momentum theorem, a substitution involving 
•other variables was found for the variables that appear in the Prandtl 
boundary-: layer equation. The substitutions reduced the boundary— layer 
equation from a partial differential equation to a nonlinear ordinary 
differential equation of the third order.. The equation was integrated 
numerically for the required boundary conditions with the aid of the 
general purpose computing system of the Bell Telephone Laboratories in 
the Langley Bell computing section. 

The reduction of the boundary— layer partial differential equation 
to an ordinary differential equation in the cases of the flat plate, the 
Jet, and the combined flat plate and Jet flow is made by finding suitable 
substitutions for the variables that appear in the Prandtl boundary— layer 
equations. By making these substitutions the velocity distributions in 
planes normal to the direction of the main flow are assumed to be similar. 
The success of the substitutions in reducing the partial differential 
equation of the boundary layer to an ordinary differential equation in 
these special eases means that the assumption of similar profiles is 
compatible with the equations describing the motion and that the similar 
profiles exist wherever the assumptions made in deriving the equations 
that describe the motion are valid. 
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The fact that the profiles are the similar profiles predicted by 
the laminar— boundary— layer theory has been verified experimentally both 
for the flow of a uniform infinite stream over a flat plate (reference 10) 
and for the flow in a round Jet (reference 11). The experimental investi- 
gation of the flow from the Jet verified the analysis (reference 9 ) for 
fll 1 regions except in the immediate neighborhood of the orifice. The 
application of the Prandtl boundary— layer equations to the present case, 
an application that implies assumptions similar to those made in analyzing 
the flow of a uniform stream over a flat plate and in analyzing the flow 
from a round orifice, is therefore expected to lead to results that are 
valid for regions not too close to the slit. 
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SYMBOLS 


nondimens ional distance parallel to surface of plate (x/a) 
distance parallel to surface of plate 

nondimens ional distance normal to surface of plate (y/a) 
distance normal to surface of plate 

nondimens ional velocity component in x— direction (uf U a ) 
velocity component in x— direction 

nondimens ional velocity component in y-direction (v/u a j 

velocity component in y— direction 

density 

coefficient of viscosity 
kinematic viscosity 

reference velocity 
reference length 

reference Reynolds number • (u a a/v) 
function of q 


function of x 
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0 


function of x 


(¥\ 


k = — ^ dT ^ 0 = A %\ 

' r 2 ~wl 

i dr) /Q 

J° 

p static pressure 

G constant 

Ci, C 2 constants of integration 

M momentum in x— direction 




reference momentum in x— direction, one side of plate ^pU a ^a^ 


P ^ 


F = 


f dt) 


i-f 


lower limit of | in 


T oo 

/dG' 
. \dfe 


•fd| 

d| / 


o-E 


a - ^V 73 
9 l/3 


7 = 


R 


1/3 


a 


surface friction f p[ ^ 

Sy 


0, 


Q- 


quantity of flow, one side of plate 
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Ujjjgjj. maximum velocity at a section 

y max distance from surface of plate to point of maximum velocity u TtlBT 

Subscripts: 

0 conditions 

00 conditions 

ANALYSIS 


at surface of plate 
at infinity 


After assuming that the flow is incompressible, that no static- 
pressure gradient exists in the direction of x-axis (fig. 1), and that 
the usual assumptions of the boundary— layer theory (reference 9) apply, 
the equations describing the motion are; The boundary-layer equation of 

motion with = 0 

bx 


s^ + f 2a = v ^ 

&C by ^2 


( 1 ) 


the equation of continuity for incompressible flow 


du bv 
bx by 


( 2 ) 


and the momentum theorem applied to the x-component of the flow 


d 


r _ 2 _ 

pu dy 


(JO 


cbc 



( 3 ), 


In order to make the equations nondimens ional, the velocities are divided 
by a reference velocity and the lengths by a reference length. The 
nondimensional velocities and lengths are 
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and' 



a 


Now equation (3) 1)60011168 


n 00 


p— p_ 

u^j/a dy 


10 


= -v 


W&u> 


a dx 


>/0 


or 


d 


P" 2. 

u dy 


JO 


dx 


1 _ 

R a 



(4) 


where 



The velocity distributions in planes perpendicular to the x— axis 
are now assumed, to be similar to one another. Velocity distributions 
are said to be similar if all the velocity distributions are given by 
one curve when u/u^ is plotted against y/ y mav . In order to make 

the assumption of similar profiles, the independent variables y and x 
are replaced by new independent variables t) and x by means of the 
following substitutions (reference 12 ): 
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where f = f(t)) and 0 = 0 (x), and 



where P = P(x). The expressions for u / u max and- for y/ y mflY for a 
fixed value of x then are 


u 




^ u=u max 


and 


l - = T 
W ( T l) u = u mflV 


All the velocity distributions in planes perpendicular to the x-axis are 
therefore given by one curve of u/u^y plotted against y/y maY and 

consequently the foregoing substitutions for y and x mean that similar 
profiles have been assumed. Whether the assumption of similar profiles 
is compatible with equations (l) to ( 3 ) which describe the motion is 
decided by whether these equations reduce to ordinary differential 
equations. 


Determination of Eelations for 0, p, and u 

The functions 0 and p are to be determined by using equations (1) 
to ( 3 ). If the values 


f 

0 


du _ 1_ df 
dy 0P dq 


and 
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are substituted in equation (4), the result is 


poo 



cbc 


i_ i_ 

R a 0 P 



or , since 


n oo 


f 2 drj 


JO 


is independent of x 


Now let 


1 dp _ 2 

(3 dx 0 dx 




L 

p 2 


than 



do 


1 dp_ 2 d 0 __E^ 
P dx 0 dx - p 2 


Equation ( 5 ) can be solved by using equation (l) to provide a relation, 
which does not involve x, between 0 and (3. After equation (l) is 
made nondimens ional, the result is 

Su du _ 1 c^U 
+ ^y R a By 2 


u 


( 6 ) 
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Equation (6) ia to "be rewritten in terms of f, <f>, (3, and q 

the relation 


it follows that 


f 

u = — 

0 


and. 


du 1 df 


du 

dx 


dy 

03 dq 

d 2 u 

_ _ 1 _ 

d 2 f • 

dy 2 

03 2 

dq 2 

U 

df dp 


V- 

dq dx 

f 


30 


0 ' 


2 dx 


From equation (2) 


where 


or 


iy py nr\ 

„ . dv , du. du Q , 

V = — dy = - — dy = - — -3 dq 

lo ^ Jo ^ Jo ^ 


dy = 3 dq 


M-p, 

t = — „ £ d „ + 

i* Jo 4,1 f J 


M 

f dq 

0 


. From 
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By partial integration 


JO 


df , 

B T~ d B 

&T\ 


= f T) - f dT| 
0 


then, the result for v is 


= $i^ + | f 


0 dx 


0. 


02 dx 0 dx 


X 'V N.P 

If values for u, v, — , — , and — - are substituted in equa- 

dx by dy 2 

tion (6), the result is 


/ 

f l _ d B dp + f_ d0 

0 \ 30 dx 02 dxy 


£lM + 

0 dx 


'0 


f d„(£-M isp> 

0 2 dx 0 dXy 


03 IB B a 0 3 2 dT) 2 


The left-band and right-hand sides of equation (7) are multiplied 
by 0 2 , and the factor 0/p 2 which then appears on the right-hand side 
of the equation is replaced by the value of 0/p 2 from equation (5). 
After further simplification, equation (7) can be written as 


1 d£ 

0 dx 

1 dp 
P dx 




d 2 f 

dr ] 2 





(8) 


The variables x and tj are independent. Consequently, a change in x 
does not change the right-hand side of equation (8) and therefore does 
not change the left— hand side either. A change in t) likewise does not 
change the left-hand side of equation (8) and thus does not change the 
right-hand side either. Both the left-hand and right-hand sides of equa- 
tion (8) are therefore equal to a constant C. Then 
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1 d0 

0 dx 
= C 

1 dp 

P dx 
or 

£0 = c dp 

0 P 

or, by integrating 

0 » cy ( 9) 

When equation (9) is substituted in equation (5), the result is 
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The expression for u contains two arbitrary constants C-^ and. C 2 . 

The constant C ± is a scale factor; the constant C 2 is a translation 

factor that determines the location of the slit along the x-axis. There- 
fore, the constant C]_ can he made unity and the constant C 2 can be 

made zero without affecting the generality of the solution in any way 
other than fixing the origin of the slit at x = 0. The expressions 
for 0, 0, and u therefore become 


and 




u = 


Kx 


2C - 1 


(2 - C) 


C 

2-C 


f < 


y 


r kx 

[_2C - 1 


(2 



The constant C can be evaluated by use of the expression for the 
momentum in the x-direction: 


or 



M 

pU a 2 a 



A00 
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or 


M 

pUa 2 ^ 



1— 2C 

2- C 




( 10 ) 


The arbitrary length a is interpreted as the distance along the 
x-axis from the origin to the point where the similar profile has been 
established. At this distance a from the origin the rate of momentum 

— O — 

flow is equated to pU & a; thus, the reference velocity U a is defined. 
When x = a or x = 1, M=M a and, consequently, equation (10) at 
x = a becomes 


1 = 



(2 - C) 


1— 2C 

2- C 



(ID 


The expression 


^ OO 

f 2 dt) 

Jo 


is a constant. 


Let 



The function f to be found mu3t satisfy this condition. Then equa- 
tion (ID becomes 


Therefore 




1— 2C 

2- C 


2K + 1 


C 


K + 2 
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Then 


K+2 

3 = x 3 (12) 

2K+1 

0 = X 3 (13) 

and 



.( 1*0 


where now 



(15) 


and 


x > 1 


It can now be shown that the requirement that C f= 2 and t ha t C / — , 

which was stipulated in determining 3 as a function of x by using 
equations (5) and (9), is satisfied. The expression that was previously 
obtained for C is 


n 2K + 1 

ITT 

From physical considerations 0 < (— \ < » and 0<R < «. There— 

V dT i Jo a 

fore, 0 < K < °°. Thus, because K > 0, 


the expression for C gives 
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Because K < °°, 


C < 2 


Therefore , 


C £ 2 

and 



In order to show that the flow appears to come from a slit^ the 
edge of the region in which the velocity u is other than zero' can he 
arbitrarily defined by choosing a fixed value of tj that corresponds to 
a small fixed value of u. Thus, let the edge of the jet region be 
given by 


1 = ^edge = Constant 


Then 


*edge 


^edge 
K+2 
„ 3 


= Constant 


or 


K+2 

y e dge = x ^ X Constant 


The flow therefore does appear to come from a slit located at 

(x = 0, y = 0). The slit can be placed at a position of x that is 

not zero by choosing a value for C 2 that is not zero. When the flat 

plate is removed, = 0 or K = 0; and then from. equation (14) 
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u 




(16) 


which is similar to the relation for u in reference 8. 


Derivation of the Differential Equation for G 

The first step in determining the differential equation for G is 
to substitute equations (12) and (13) in equation (7). The equation of 
motion, therefore, becomes 


2K ± l f 2 
3 


+ 





1_ dff 
E a dq 2 


( 17 ) 


The assumption of similar profiles has reduced the partial differential 
equation of the boundary layer to an ordinary differential equation; 
thus, an indication is given that the assumption of similar profiles is 
consistent with equations (l) to (3) which describe the motion. Now. 
let 


1TJ 

f dt|=J 

Jo 

where F = F(q). Then, because 



df = d^F 
dT l ” dq 2 


and 

d 2 f = d^F 
dq 2 dq^ 
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equation (17) becomes 


1 _ d^F 1 - K dfF 2K + 1 /dF \ 2 

E a d ,3 + 3 dT1 2 + 3 {^) = 


Now, let 



and 


where G = G(|). Thus, 


a 



F = -fi 


and 


dF ay dG 
drj _ 3 dl 

d 2 F = a^2i dfG 

d q 2 9 d | 2 



g3y d3G 
27 d|3 


Equation (l8) then becomes 


(18) 


^2 + (l_K)G'S!= + (l + 2K) (mf . 0 

d|3 d| 2 '«/ 


(19) 
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where 

g g y f d g t>\ _ a 3 / d g G\ 

9 E a\d|2/ 0 9B a g Vd£ g / 0 

Now, let 



a = 9 1 / 3 E a 2 / 3 


Then, equation (19) becomes 



The conditions to he satisfied by G 


are, at 1=0, 


G = 0 



and, at £ = oo. 


(20) 


G = Constant 


Integration of the Differential Equation for G 


The differential equation for G, equation (20), is integrated 
numerically by use of the method of reference 13. The numerical inte- 


gration was started by computing the values of G, 


dG 



and 


for the first five intervals of | by means of a Taylor's series.. 
Taylor's series that is valid for small values of I is 


d| 3 

The 
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From tli© nature of the problem each term of the equation is continuous 
in th© interval from 0 to the equation may therefore be integrated 

directly term by term. Then 




= G 


d%r_ r 

i^Jo Jo 


d% dG 
d| 2 'll 



d^G dG . t 
— r — d £ 
d| 2 d l 


Integrating 


n°<> o 
dfG dG 

Jo M 2is 


d| by parts gives 



a% as 
as 2 « 


3"i 00 

'«'J 0 


1° 

Jo 


“i^as 

at 2 



d 2 G dG 

d| 2 d i 


dl 


Therefore 


p°o o 

d G dG 

Jo 


as 


0 



d| = 0 


and thus 
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Now, integrating 


G 2 d% ^ by parts gives 
d| 2 


n oo 


g 2 d? = o 2 22 

dt 2 d£ 

0 d 5 


_,oo n oo 


0 


— 2G — d| 
d£ d| 


o°° 


Jo 


2 O(^) « 


Then equation (21) becomes 


-2 


1 - 


Vdl 2 


noo 


JO 


G ( — 
,d| 


d| + 


1 + 2 


d%\ 

^ 2 /o 


G^f 

\di y 


d| = 0 


or 


v d ^ 2 /0 


- 1 


noo 


Jo 


\2 

g f^y di = o 


( 22 ) 


The quantity G is always positive, except on the plate where it 

'd 2 ^ 1 


is zero; therefore, equation (22) cannot be true unless 


M 2 )o ' k ' 


Consequent ly, 


/ a 2 G\ 1 


^ 2 Jo ^ 

the numerical integration of equation (20). 


is taken as the starting value of 


dgG 

di 2 


in 


The equation was integrated by use of the "fourth approximation" of 
reference 13 for an increment in |, A| = 0.10; the numerical integration 
began at % = 0.50. In order to determine whether the increment in 
£, Ag =0.10, was sufficiently small, the equation was also integrated 
for an increment in g, Ag = 0.05; the numerical integration began 
at | = 0.50. Some of the results from the integration with Ag = 0.05, 
the more accurate integration, are given in table I. A comparison of 
some of the values of G and its derivatives obtained by integrating 
with Ag = 0.05 and with Ag = 0.10 is given in table II; the values 
at g = 0.60 computed by the Taylor’s series are also included for 
comparison with those computed by Falkner's method (reference 13).. 
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The results of the two integrations agree to five significant 

figures up to £ = 8.50 at which value i s less than 1 percent of 

d| 

its maximum value. From | = 8.50 to | = I3.IO, at which value the 
computation for A| = 0.10 was stopped, the agreement between the 
computations for the two different increments in | becomes poorer. 

The values of G, however, are the same to five significant figures for 
both computations. At | = 13.00, G has already attained its asymptotic 
value, a value that is necessary for the computation of v at | = 00. 

(See table I.) 

As a further check on the accuracy of the numerical integration, 
the following condition was used: 


T 00 

f% = 1 

Jo 


which is the same as 


noo 


do 


:)% . i 


( 23 ) 


The values of 


PI 2 

/dG\ 

\ d| obtained by using Weddle's rule (reference 14) 


do 


with A| = 0.05 are given in table III. It is seen that the foregoing 
condition (equation (23)) is confirmed to five significant figures. In 

poo 

order to determine whether the integral 
with sufficient accuracy by the integral 


do 


( — ) d| is approximated 
\ d| / 

the assumption 


112. 30 

1 


\ 

JO 



was made that for large values of | the curve of — plotted against | 

d| 

could be approximated by the function 
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The constants A and B were evaluated In the range of | from 1 = 9. 
to | = 14.95* The const ant B was observed to increase slowly. For 


dG . -B| 
— = Ae 
d| 


the integral from 1-^ to <» is 



= (o-oooo^^Tio) 2 

2(0.9452) 

= 0.4425 X 10 -8 


where 


i L = 12.30 


and 


= 0.00009145710 


B = 0.9452 


dG 


Because — decreases as increases and because B increases 


\ d| Is given accurately to 

r\ \ / 

P12.30 


Jo 


five significant figures by the integral 


do 
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T oo 

In order to determine whether the integral (— \ d| is approx i— 


0 

14.95 


mated with sufficient accuracy. by the integral 




the same 


dG 


function for — for large values of £ was assumed; that is, for 


it 


the integral from to <» i s 


where 


and 



_ o. 000007433791 

0.950682 
= O.OOOOO782 

t L = 14.95 


dG 


O.OOOOO743379I 


B = 0.950682 



26 


MCA TN No. 1644 


_ d.G 3 

Because — decreases as increases and because B increases 

n oo 


as ^ increases, adding 

Jl4.95 

does not change its value to five significant figures. 


— d| to the value of G at | = 14.95 


n °° 


The last value in table III for 


JO 


^dG 


dGV 

<4/ 


d| and the asymptotic 


value of G are therefore both unchanged* by extension of the curves 
to i = 00. The results of the computation for A| = 0.05 are conse- 
quently taken as correct to five significant figures for 0 ^ | £ 8.50 

and the asymptotic value of G is taken as 1.2599. The curve of — 

d| 

plotted against I is shown in figure 2. 


Determination of* Final Expressions 

Nov that the values of G, and have been tabulated, 

^ d| 2 d|3 

the completion of_the solution of the problem requires that expressions 
for u, u, v, v, Q^, t q , m, and | be determined. 


In order to obtain the expressions for u and u, use 


dF 

u = £ = -*]_ 
0 2E+1 

„ 3 


0J2; dG dG 

3 4 _ 9 2 /^ 1/3 dg 
jl/2 ' 3 1 *1/2 


(24) 


The expression for B a is 


R a = 


V a 


1 


3a 

pa 


\l 


U£ 

pv 2 


therefore, u may now be written as 


^ \PV 2 / X 1 / 2 


(25) 
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The expression for u is 



-■ //- _ x / 2 dG 

1 « 

3 VpvV ^ \Ip S 


or 


u 


til M 1/3 
3 w 


- 1/6 

8. dG- 

U7 2 d| 

x 


(26) 


The u component of the velocity therefore varies Inv ersely as the 
1/2 power of the distance from the slit. 


The expression for v is 


= £i4S + I f dr / £_ M _ ± 


0 dx 


JO 




;2 dx f dx 


where 


f = ££ 

3 d| 

7 

^ , K+2' 

x 3 

K+2 

3 = x 3 


2K+1 
0 = x 3 
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Pn 


f dr] = yG 


I/O 


and 


I = 


at) 

3 


Then 


v = -L- f *+-£ dG K-l > 

Kf2 A 3 d| 3 j 

x 3 


and, after the appropriate substitutions for 7 and K have been made. 


_ 9 1 / 3 1 A , dG _ G n 

"b^/s^A^ « 


or 


v = 


_ 9 i/3 


(»«-) 


.PV / 


(27) 


The expression for v is 


v = vU Q = 


9 1 /3 1 . 

k -3A 4 5 


- 3/u / 3 S iS_ G ^|^ 

pa 


■PV 


2 J 


or 


X 


(28) 
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The magnitude of the vertical velocity v thus varies inversely as the 
3/4 power of the distance from the slit. Whether the vertical velocity 

is away from or toward the plate depends on the sign of the term 3£ — _ G 

in equation (28). 


When | — » °°, it can he shown that | — — >0 because G — > Constant 

d| 

as | — >oo. Therefore, 


or 


9 V 3 / V \ 1/3 IT 1712 „ 
**■ \ p / _3A ‘ 


- = _ (1.2 999 )9 1 / 3 /MaVX 1 / 3 i 1 / 12 

4 -UJ 57T 

z 


(29) 


The flow from the slit thus induces a velocity towards the plate at 
sufficiently large distances above or below the plate. 

, , T 11 ® q - uanti ty of fluid passing through an imaginary plane normal to 

the plate and extending to. infinity in one direction is given as 


noo 


Qoo = u dy = U a i 

Jo 


u dy 


1)0 


and, after the appropriate substitutions for U & , u, and y have been 
made. 


or, finally. 


q m = 9 1/3 x 1/4 ^ J /3 v 1/3 j ^ 


d| 



( 30 ) 
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where 


&00 = 1.2599 


The quantity of fluid, passing through a plane normal to the plate per 
unit time thus varies directly as the 1/4 power of the distance from 
the slit to the normal plane. 

The expression for the surface friction t q is 



^ U a 



1_ /df\ 

a 0P \drJ o 

^a. o?z /d%\ 


*°a 5 ?A k 
l _5/4 a 



and, after the appropriate substitutions for U , R_ , and 

ct- 7 8. " 

have been made. 



T 


O 



0 
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or 


Ma a lA 1 

T ° ' %p J>A h 


1/4 

Ma a 
4 _5/4 

X 


(31) 


The surface friction is thus independent of the viscosity of the fluid 
and varies inversely as the 5/4 power of the distance from the slit. 

The shearing stress is independent of the viscosity of the fluid 
not only at the surface hut also between any two adjacent layers in the 
flow. For the flow from a slit without the presence of the plate (refer- 
ence o), the shearing stress between any two adjacent layers of fluid is 
also noted to be independent of the viscosity of the fluid. The shearing 
stresses between adjacent layers are independent of the viscosity in both 
the problem presented herein and in the flow from the slit without the 

plate, because in both cases the velocity derivative — is inversely 

Sy 

proportional to the viscosity 4 of the fluid. 

The ratio of the rate of flow of momentum parallel to the plate at 
distance x from the slit to the rate of momentum flow at the reference 
distance a from the slit can be obtained from equation ( 10 ) as 


M_ 


M 


a 1 / 4 


fT 2 - 1/4 

P^a a x 


(32) 


The rate at which momentum parallel to the plate passes through a plane 
normal to the plate therefore varies inversely as the 1/4 power of the 
distance from the slit to the normal plane. 

The expression for | can be obtained as follows: 


or 


= OTL = 9V 3 Ra 2 / 3 y = 91 / 3 /gM 1 / 3 


3A 




r 3A 


I = 


9 1 / 3 /aMaX 1 / 3 y 1 

3 \v 2 pj 3, 3 /V A 


= 9 1 / 3 /Ma \ 1/3 ^L /12 _y_ 




JA 

X 


(33) 
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The width of the disturbed region, given by y for | = i edge » thus 
increases as the 'i/k power of the distance from the slit. 

Validity of Solution 

A restriction on the region in which the solution is valid can be 
obtained by noting that for the boundary— layer equations to be valid the 

value of ^ must be small. By use of equations (26) and (28) the ratio 
u 

becomes 


. £ (»;-°) 


J/6 


3 V 2 


p v 


dfi 

-1/2 d| 

X 


or 


V = 3 f pvf 


vV3 


x a* 



Ma y 

Large, values of give small values of — , and so the region in 

PV 2 u 

Ma 

which the solution is valid is enlarged when is increased. 

PV 2 

By use of the complete equations of motion instead of the boundary- 
layer equations the difference between the static pressure on the plate 
and the static pressure very far above the plate is given approximately 
at a given value of x by 


/ Ap \ 

.. / o\2/3| 

H 

O 

H 

1 ? 

1 r 

- 2 

\P u max /_ 



l!> I 13 
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and the pressure gradient in the direction of x is given approximately 

t>y 


_5 _1_ 

_l/6 _3/2 
a x 




V/hen it is noted that v is a small quantity, the approximation 
that the static pressure is constant throughout the entire field which 
was made in equations (l) and (3) is Justified "by the solution.. All the 
terms neglected in the boundary— layer approximations can likewise be 
shown to be truly negligable on the basis of the obtained solution. All 

the approximations become more accurate as x increases and as - a . 

pv 2 

increases . 


because — J where aj_ is an arbitrary value of a and 


The constant a seems to appear in a rather arbitrary manner but, 
Ma ^ai\ 1 A 
Mai = 

is the momentum corresponding to the distance 5^ from the slit, it can 
be shown that a change in a does not change u, v, or t q pro- 

vided that x and y are kept fixed. 

When the plate is removed, the flow becomes the same as the flow 
in a two-dimensional laminar -Jet : . This problem has been treated in 
references 7 and 8. The removal of the plate in the present problem 

/ d 2A 

makes ( | = K = 0. Equation (20) then becomes 

W Jo 


d% „ d 2 G /An ' 2 
+ G + 

d| 3 d| 2 


(**f. 0 

Ws ) 


When the conditions that 


G = 0 

d^ 


= 0 


d| £ 
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at |=0 and that 


"loo 

Jo 




are used, the solution, which in this case is obtained in closed form, 
is 


G = tanh — 
2 


(34) 


Because the momentum Ma is the momentum for one side of the plate in 
contrast to the momentum used by Bickley Mg (reference 8) which was 
the total momentum flow from the slit and therefore equal to twice Mg_, 
the result £ = 2|g is obtained. When this relation between the value 
of | of the present analysis and the value of |g in reference 8 is 

used together with equation (34), exactly the same relations are obtained 
for all the- quantities of interest as were obtained by Bickley 
( reference 8) . 


CONCLUSIONS 


The laminar flow of a slightly viscous incompressible fluid that 
issues from a slit and passes over a flat plate is investigated in a 
region far enough from the slit for the boundary— layer equations to be 
valid. The results obtained are that the velocities parallel to the 
plate vary inversely as the 1/2 power of the distance from the slit, 
that the velocities perpendicular to the plate vary inversely as the 
3/4 power of the distance from the slit, that the width of the disturbed 
region increases as the 3/4 power of the distance from the slit, that 
the surface friction is independent of the viscosity of the fluid and 
varies inversely as the 5/4 power of the distance from the slit, that 
the rate at which momentum parallel to the plate passes through a plane 
normal to the plate varies inversely as the l/k power of the distance 
from the slit to the normal plane, and that the quantity of fluid 
passing through a plane normal to the plate per unit time varies directly 
as the 1/4 power of the distance from the slit to the normal plane. 

Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
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TABLE III 


PS 


VALUES OF INTEGRAL 



OBTAINED BY USE 


OF WEDDLE'S RULE FOR A| = 0.05 


£ 

l 

pi 

\d i 

0 

5 

l 

pi 

(£)% 

0 

0 

0 

6.30 

0.3329781 

.30 

.0005619072 

6.60 


.3331301 

.60 

.004462288 

6.90 


.3332176 

.90 

.01476475 

7.20 


.3332678 

1.20 

.03369614 

7.50 


.3332965 

1.50 

.06194261 

7.80 


.3333129 

1.80 

.09818613 

8.10 


.3333222 

2.10 

.1392639 

8.40 


.3333275 

2.40 

.1810349 

8.70 


.3333305 

2.70 

.2195948 

9.00 


.3333322 

3.00 

.2522585 

9.30 


.3333332 

3.30 

.2779223 

9.60 


.3333338 

3.60 

.2968261 

9.90 


.3333341 

3.90 

.3100159 

10.20 


.3333343 

4.20 

.3188173 

10.50 


.3333344 

4.50 

.3244822 

10. 80' 


.3333345 

4.80 

.3280249 

11.10 


.3333345 

5.10 

.3301909 

11.40 


.3333345 

5.40 

.3314920 

11.70 


.3333345 

5.70 

.3322630 

12.00 


.3333345 

6.00 

.3327151 

12.30 


.3333345 
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Figure 1.- Flow from a slit at 0 into the XY -plane with a plate along OX. 





